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The research described in this report was prepared by Honeywell Inc. , 
Minneapolis, Minnesota 55413, under NASA Langley Research Center Con- 
tract NASI -154 86. This work was directed by the Flight Mechanics and Con- 
trol Division of the NASA Langley Research Center and was administered by 
Mr. Henry J. Dunn of the Flight Mechanics and Control Division. Special 
thanks to Mr. Jarrel Elliot and Mr. Irving Abel for their continued support 
toward this contract. 

The technical work reported in this volume was conducted by the Research 
Department at the Systems and Research Center of Honeywell Inc. Dr. A. F. 
Konar was the Honeywell Program Manager. Mr. C. R. Stone and Dr. J. K. 
Mahesh were the principal investigators on this contract. They were assisted 
by Mr. P. D. Hausman and Dr. W. L. Garrard. This report covers the work 
done from August 1978 to August 1979. 

The work under this contract was reported in three volumes entitled 
"Active Flutter Control for Flexible Vehicles. " 

Volume I Final Report 

Volume II Appendices 

Volume III KONPACT II Program Documentation 
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SECTION I 
INTRODUCTION 


The general objective of this study is to develop methodology for rapid 
design of active control systems for flexible vehicles. Active control can 
help meet the goals of weight reduction and improved efficiency demanded of 
future flight vehicles. The methodology for active control synthesis is based 
on linear quadratic gaussian (LQG) theory. This is a viable method for 
design of multi- input multi-output controllers to meet multiple, often con- 
flicting, performance requirements. This methodology also gives the active 
control designer the capability to influence the geometry of the vehicle to 
realize the benefits of control configured vehicle (CCV) design. 

The scope of this study is to develop an active flutter control methodology, 
including surface/sensor placement, that will apply to a wide class of flexible 
vehicles and demonstrate it by designing a flutter controller for DAST ARWl 
wing to increase the flutter speed by 20 percent. This study also included 
the development of algorithms and computer programs for flutter modeling 
and active control design procedures. These have been integrated into NASA 
owned KONPACT software (ref. 1). The resulting software is called 
"KONPACT II-Computer Programs for Active Control Technology. " 

The work under this contract was reported in three volumes entitled 
"Active Flutter Control for Flexible Vehicles. " 

Volume I. Final Report 

Volume II. Appendices 

Volume III. KONPACT II Program Documentation 

Volume I reports tne active flutter control methodology and the summary 
of the technical results obtained under this contract. 


An overview of the Honeywell's approach to the synthesis of active 
flutter control for flexible vehicles is given in Section II. A brief descrip- 
tion of KONPACT II software is also described in this section. 

In Section III, the dynamic modeling of the DAST ARWl wing for optimal 
flutter control design is presented. This is an automated modeling procedure 
which minimizes human error due to data handling. 

The results of optimal state feedback synthesis for the DAST ARWl wing 
are presented in Section IV. The control surface optimization results are 
also discussed in this section. 

Section V presents the robust Kalman estimator design procedure and its 
application to the flutter control of the DAST ARWl wing. Here the robust- 
ness is characterized as gain and phase margins of the overall system. Also, 
the sensor optimization/selection methods and results are discussed in this 
section. 

Practicalization, of the optimal flutter controller with robust Kalman 
estimator, is presented in Section VI. The evaluation of the final flutter 
controller at different mach numbers and dynamic pressures is also presented. 

Section VII contains the conclusions and recommendations. 

The detailed technical results during the course of this study were docu- 
mented as Customer Engineering Letters. These are collected in the form 
of appendices and presented in Volume II. 

, Volume III contains the description of the computer programs developed 
during this study. The user's information is also provided in that volume. 
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SECTION II 

ACTIVE FLUTTER CONTROL SYNTHESIS METHODOLOGY 


This section presents an overview of Honeywell's approach to synthesize 
active flutter controller for flexible vehicles, The main design objective is 
to increase the flutter free air speed by at least 20 percent (44 percent in 
dynamic pressure) at Mach number 0. 9. 

Honeywell's design procedure using LQG theory consists of three parts: 

• Optimal state feedback synthesis 

• Robust Kalman estimator synthesis 

• Practicalization of the resulting optimal controller. 

First, the flutter design model, for the DAST ARW1 wing, is described. 
This is followed by the step-by-step description of the design procedure to 
synthesize a practical flutter controller. A brief description of the soft- 
ware tool KONPACT II is given at the end of the section. 


The Flutter Design Model 

The flutter system consists of a half wing, an aileron, actuator, and 
a set of sensors. The basic modeling task is to obtain a state space repre- 
sentation for the flexible wing with unsteady aerodynamics (ref. 2). The 
modeling procedure is discussed in Section III. A block diagram of the 
overall flutter system is shown in Figure 1. The white noise input p , is 

o 

shaped by the gust filter dynamics to produce the effect of gust (wind) force 

on the wing. The control input, u c , moves the control surface through the 

actuator dynamics to coritrol the fluttering of the wing. The output of the 

sensors is m , 
s 
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Optimal State Feedback Design 

In this control synthesis step, it is assumed that all the states (x) of 
the flutter design model are available. The control law is obtained as a 
set of feedback gains. The optimal state feedback control law represents 
the best that is achievable. It is guaranteed by theory to have good robust- 
ness properties (ref. 3). As a more of the practical requirements are en- 
forced, the performance of the overall system deteriorates. So the optimal 
state feedback control should be synthesized with stringent performance 
requirements. The details of this design procedure are the subject of 
Section IV. 


Robust Kalman Estimator Design 

This step eliminates the need to have all the states of the flutter design 
model. It is assumed that enough sensors are placed on the wing so that 
the states of the flutter design model can be estimated. The flutter control 
law is then obtained by multiplying the estimated states by optimal state 
feedback gains (obtained in the previous step). Here the designer has the 
freedom to choose the number and location of sensors. Also, there is free- 
dom while designing the robust Kalman estimator to improve the robustness 
property of the overall system at the expense of rms performance. These 
issues are discussed in Section V. 

Practicalization and Performance Analysis 

The robust Kalman estimator together with the optimal state feedback 
gains constitute a flutter controller that is implementable. However, the 
cost and effort of implementation can be significantly improved by approxi- 
mating the above flutter controller by a low-order one. The two procedures 
suitable for this task are residualization and frequency response matching 
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f 


(ref. 4 and 5). Once the flutter controller is practicalized, it should be 
evaluated with a detailed performance analysis at various mach numbers 
and dynamic pressures to assure its operation throughout the flight enve- 
lope. The details Of this design step are discussed in Section VI. 

KONPACT II Software for Active Control Design 

KONPACT software (ref. 1) was developed to automate the modeling 
and design steps of active control to minimize the error of data handling. 
The software uses advanced computational techniques to perform system 
modeling, modern control synthesis, analysis and design of active control 
systems. This software was updated with respect to modeling procedures 
and design methods developed during this study. The resulting software 
is called KONPACT II. Figure 2 shows the program organization, for 
more details on the software the reader is referred to Volume III. 
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Figure 2. - Program organization of KONPACT II. 













SECT ION\ III 

dynamic modeling for flutter control design 

The flutter system consists of a half wing, an aileron, actuator, and 
a set of sensors. 

The wing is a dynamically scaled representation of a transport-type 
research wing. It is scaled to flutter within the operational limits of the 
NASA Langley transonic dynamic tunnel (TDT). It is provided with a hy- 
draulically actuated trailing edge control surface. Figure 3 is the model 
geometry of the wing. It has an equivalent aspect ratio of 6. 38 and a leading 
edge sweep of 44. 32 deg. The trailing edge control surface is 20 percent 
of the local wing chord and is located between the 76. 3 percent and 89. 3 per* 
cent semispan stations. 



Figure 3.- Model geometry of the DAST ARWl wing. 
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Modeling the Dynamics of the Wing and the Aileron 

Using analysis technique described in (ref. 6) the equations of motion 
can be written as 



( 1 ) 


where 

= Vector of generalized coordinates for flexible modes 
= Vector of generalized coordinates for control surfaces 
= Vector of generalized coordinates for gust inputs 
M = Generalized mass matrix 
K = Generalized stiffness matrix 


Q(k) is the unsteady aerodynamic influence coefficient matrix evaluated 
at various reduced frequencies k given by the expression 


k 


DC 

2 V 


( 2 ) 


where 

to = Circular frequency 
c = Reference chord length 
V = Velocity of the vehicle 


and q is the dynamic pressure and is given by the expression 

q = 1/2 p V 2 (3) 

where 

p = Density of air 
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For this study, NASA, LaRC provided a NAS.TRAN finite element model of 
the wing. The first 10 elastic modes, generalized masses, and frequencies 
were provided. The modes cover a frequency range of 5. 23 to 118. 15 Hz. 
The unsteady aerodynamic influence coefficient matrix data was provided for 
eight reduced frequencies ranging from 0 to 1. 2. The mode shape data at the 
various nodes of the finite element model was also provided to compute the 
sensor (accelerometer) output equations. 


The state space flutter model is obtained by the least square approxima- 
tion of the unsteady aerodynamic forces. This approximation has been shown 
to 1 give accurate models for flutter analysis and design (ref. 6 and 7). The 
unsteady aerodynamic force matrix Q(s) is approximated by a rational poly- 
nomial, in the Laplace transform variable s, given by 


Q(s) 





+ A, 



L 


+ E 

m = l 


D s 
m 



(4) 


where Aq, A^, A g, D , m = 1,L are the aerodynamic coefficient matrices 
and k , m = 1, L are the aerodynamic lags. The original form of this 
approximation was suggested by R. T. Jones (ref. 2). 


The values of are selected, within the range of reduced frequencies 
over which Q(k) has been calculated, to minimize the error of approxima- 

A 

tion. Usually only a few values of k m near the reduced flutter frequency are 
required. The aerodynamic coefficient matrices are computed by least- 
square fit of the aerodynamic force at k., i = 1, N. The least-square proce- 
dure and the conversion of the resulting s -plane approximation to state space 
description by using minimal realization technique is described in Appendix 
A (Volume II), Figure 4 is a block diagram of the state space realization 
of the flutter model. The root locus at mach number = 0. 9 as the dynamic 
pressure is varied is given in Figure 5 (the flutter model has L = 4 and k. 

A A A X 

= 0. 2, kg = 0. 4, kg = 0. 6 and k^ = 0. 8). Note that as the dynamic pressure 
is varied the frequencies of modes 1 and 2 tend to coalesce exhibiting a 
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Figure 4. - State space realization of the flutter model. 
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typical bend ing - torsion type flutter behavior. The results are in agree- 
ment with the experimental flutter studies conducted by NASA in the Langley 
TDT (ref. 8). 

The Actuator Dynamics 

The trailing-edge control surface is driven by an electrohydraulic servo 
actuator system. If the aeroelastic effects and the mass of the aileron are 
neglected, the deflection of the actuator will be the same as that of the aileron. 
NASA measured the actuator transfer function experimentally during wind 
tunnel tests (ref. 9). The following transfer function approximates the mea- 
sured transfer function: 



where u is the control signal and 6 is the control surface deflection. The 
aileron rate and acceleration are obtained by implementing the actuator 
dynamics given by equation (4) as shown in Figure 6. 

Based on the maximum control surface displacement and rate, NASA 
specified the maximum RMS deflection and rate to be 4 deg and 250 deg/sec, 
respectively, for an RMS gust velocity of 0. 3048 m/s. In addition, a control 
system frequency response attenuation boundary was specified to meet the 
actuator bandwidth requirement. Figure 7. 

The Gust Model 

The RMS control surface displacement and rate are determined by using 
a second- order Dryden wind model. 
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The Von Karman power spectrum to describe the vertical component of 
one-dimensional atmospheric turbulence, (ref. 10) is 


a 2 L 

0 (y) = 

ug rrV 


where 

L = 30. 48 m 

V = 137, 16 m/s 

a = 0. 3048 m /s 

wg 


1 +' 


1. 339 


1 + ■ 


1. 339 


Lu 

V, 


Lto 
Vi , 

oii/i6 


( 6 ) 


The second -order Dryden wind model to represent the Von Karman 
power spectrum (ref. 10) is 


u ■■ o 
g = 

g 


\fei i+ \F¥ 


(7) 


1 +is 


Substituting for L, V and 


u 


g - 2, 07 (s + 2. 6) 


n 


g 


(s + 4. 5 


( 8 ) 


A first-order Dryden model is obtained by approximating equation (8) 
over the frequency range from 0 to 100 rad /sec. It is given by the transfer 
function 


W g = 2,5926 

~ (s + 3. 3607) 
g 


(9) 


The two wind models over the frequency range of interest are compared 
in Figure 8. 
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Figure 8. ~ Comparison of two dryden wind models. 
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The Design Model 


The design model is obtained by combining the wing, the actuator and the 
gust model and defining the design responses for LQG design (ref, 4), 

Figure 9 is a block diagram of the design model. The requirement on the 
design model is that it should be of as small an order as possible so that the 
design costs are reduced. Residualization is a reduction procedure well- 
suited to approximate the high-frequency modes (ref. 4), In the case of the 
flutter model, the parameter L, which denotes the number of aerodynamic 
lag terms, provides another natural procedure for reduction of order of the 
design model. In this study, a reduced-order design model was obtained by 
retaining modes 1, 2, 4, 5 and 6 and using L = 1, The order of this design 
model is 20 states. This reduced-order model preserves the open-loop 
flutter characteristic to within 5 percent. For the final performance analysis, 
a model was obtained using modes 1, 2, 4, 5, 6, 7, 9 and 10 and using L = 4. 
The order of this model is 53 states. The results presented in Section VI 
show that the performance values are within 15 percent. The reduction in the 
design cost is quite significant. 







SECTION IV 

OPTIMAL STATE FEEDBACK DESIGN AND 
CONTROL SURFACE OPTIMIZATION 


Optimal state feedback synthesis is the first step in the active control 
design methodology to obtain practical control laws. Three sets of results 
are presented for this step in this section. They are: 

• State feedback design 

• Control surface optimization, and 

• Insensitive design 

In addition, the results of control-observe-ability are also presented. 

The state feedback design results with nominal control surface are 
presented first. They show the best performance achievable at the design 
condition. Both first-order and second-order Dryden wind models are used. 
The second-order wind model enforces response values that are 40 percent 
larger. Fifteen controllers (two with a first-order wind) are synthesized, 
all controllers meet performance requirements, one with each wind model 
is chosen for subsequent design efforts. 

The results of two control surface optimization problems are presented 
next. For the first, the aileron position is optimized over those with a 
semi-span of 0. 130 (the same as the nominal). For the second, the optimi- 
zation is for any semi-span within the allowable range of 0. 660 and 0. 937 
of semi-span. For the first case, the best aileron position is as far out- 
board as possible. Reduction of RMS values of aileron deflection and rate 
of 3. 21 and 1. 73 percent are achieved. For the second case, the best 
aileron is that for maximum semi-span. Reductions of RMS values of 
aileron deflection and rate of 45. 93 and 48. 91 percent are achieved. 
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Ins ehs Rive-design results are then presented. The uncertainty weight- 
ing procedure is used to develop constant controllers that provide good 
control in spite of large variations in model coefficients that are a function 
of dynamic pressure. The price paid is a small increase in RMS response 
activity. 

Finally, an assessment procedure for the quality of controllability and 
observability is discussed. It is better than the rank test. The quality of 
controllability is reasonably well graded. The theory is also applied to 
accelerometer placement and is found to be of value. 


State- Feedback Design 


The performance objective is to stabilize the wing in the flutter condi- 
tion at 0. 9 Mach and a dynamic pressure of 7. 66 k Pa. The final control- 
ler should have RMS response values for aileron (6) and aileron rate (fi) 
less than 4 deg and 250 deg/s for an RMS wind of 0. 3048 m/ s. At the full- 
state feedback design step the control activity should be about 35 percent 
smaller to allow for the performance deterioration at the subsequent design 
steps. 

The plant model is: 

x = Ax + BjU + Bgri ^®) 

r = Cx + Du (11) 

where 

x. = n i (i = 1, 5) flexure mode deflections in m 

x. = ri. c (i = 6, 10) flexure mode rates in m/s 
l 'l- 5 

(i = 11, 15) unsteady aerodynamic states 
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x i (i = 16, 18) actuator states 
x^g = 6 Aileron deflection in deg 
x^ 7 = fi Aileron rate in deg/s 

X last = g ust velocity (scaled for 0. 3048 m/s) 

u = Scalar control input 
rv- Unit white noise 


The model is of 19-order or 20-order depending on whether the first-order 
or second-order (Dryden) wind model is used. 


The full-state feedback control law 


u = Kx 


( 12 ) 


is determined from 

j = Min E f r . Q r *j = E j x > Qx + pu 2 } (13) 

We present two sets of results with the first order wind model: Sirisena 
method and baseline. 

With the Sirisena (ref. 11) a response transfer function r/u is construct- 
ed by deleting the wind and requiring the real parts of the zeros to be -1. 01* 
(open-loop poles; except the first-order actuator pole). 

Then 

r/u = c' x/u 
Q = cc' 


(14) 

(15) 
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As p -* o, the closed-loop poles approach the numerator zeros, and 
the actuator pole goes to the wind state remains invariant. Experience 
shows the desired asymptotic properties of this "minimally acceptable" 
feedback controller are achieved at nonzero values of p with acceptable 
values of pole position. Volume II has Q and p values. 

Table 1, column case 0, lists the nonzero diagonal elements of the weigh 
ing matrix Q for the baseline case. Columns S and 0 show that both the 
Sirisena and baseline cases have highly acceptable RMS responses with the 
Sirisena being slightly better. Actuator pole positions are essentially at 
their open-loop value; the Sirisena real pole is slightly stiffer. Both are 
viable candidates for further development but we will subsequently (Sections 
V and VI) use only the baseline. 

Figure 10 shows that both the Baseline and Sirisena controllers pro- 
vide effective stabilization over a dynamic pressure range from 3. 83 to 
8. 62 k Pa (at MachO. 9). 

Columns 1-13 of Table 1 present results with the second-order, Dryden 
wind model. Cases 0 and 2 use the same quadratic weights; the RMS re- 
sponses are 40 percent greater with the second-order wind model. 

For case 3 there is weight only on control. Kwakernack and Sivan 
(ref. 12) show this is also a "minimally acceptable" optimal control; it 
reflects the unstable poles. 

For cases 1-5 there are no weights on 5 or 5. Actuator poles are near 
their open-loop values but RMS response values (particularly for aileron 
rate) are too close to their limit values. We therefore weighted 6 and/or 
6 in cases 6-13 to reduce RMS control responses. We select case 8 as a 
baseline for further design effort. 


22 


































8.618 7.661 



Closed loop designs at q = 7.661 




As is usual in the application of LQG methodology, we find here that 
its ease of use in providing tradeoff data to justify its use. In the present 
synthesis process, it is also the first step in the synthesis of practical 
controllers. 


Control Surface Optimization 

Given a scalar y of control surface positions, we wish to determine 
the optimal value of y and the feedback gain matrix K. 


u = Kx (16) 

to perform 

j _ Min Min E {r , Qr1 ( 17 ) 

y 

and Q is given by the case 8 values of Table 1. We use a modification of 
the Johnson algorithm (ref. 13) as is discussed in Appendix P of Volume II. 

To determine the optimal control surface position we start with an 
initial y and successively perturb y to get reductions in J until a minimum 
is achieved. Because of strong monotonicity, y is simply incremented to 
the boundary; we do compute the one-sided difference approximation to the 
gradient of equation (17) to be assured that we do not skip over a local 
minimum. 

Table 2 summarizes the results. The first row entry is with the 
nominal aileron and is the same as case 8 of Table 1. Row 3 is with the 
same aileron position but employs a polynominal fit model for the aileron 
aerodynamic data. There are small unimportant differences. 

Rows 2-5 establish that for ailerons whose span is 0. 130 of the wing 
semispan that the most outboard allowable position is best. Our first itera- 
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TABLE 2. - CONTROL SURFACE POSITION OPTIMIZATION 


Aileron position 

Part of semispan 

Inboard 

Outboard 

0. 763 

0. 893 

0. 753 

0. 883 

0. 763 

0. 893 

0.783 

0. 913 

0. 807 

0. 937 

0. 730 

0. 937 

0. 660 

0. 937 


Poly- 

lomina 

fit 


RMS 


8 


Deg Deg/sec 


2.579 189.4 

2.620 192.0 

2. 594 190. 5 

2. 584 188. 5 

2.510 187.2 

1.717 124.2 

1.428 97.7 


2. 638 0. 5529 x 10 

2.671 0. 5668 x 10 6 

2. 644 0. 5576 x 10 6 

2. 599 0. 5438 x 10 6 

2. 561 0. 5338 x 10 6 

1.756 0. 2981 x 10 6 

1.462 0. 2308 x 10 £ 


0. 9047 x 10^| 3 

0. 1434 xlO 5 3 
0. 2582 xlO 6 3 
0. 5753 x 10 
0. 8458 x 10 


2 2 
3 

4 4 

5 5 

6 6 





















tion (row 2) was to move the aileron inboard, 6, 6, u, and J all increased. 

The one-sided difference is an excellent approximation to the differences in 
the values of J. This assures us of monotonicity over inboard shifts in 
aileron by increments of less than 0. 01 semispan and assures us reasonably 
well that inboard shifts are undesirable. 

We then moved the aileron outboard from its nominal position in two 

* 

steps (rows 4 and 5) to the maximum allowable. The values of 8, 6, u, and 
J are monotone decreasing and DJ's are again a good approximation to the 
changes in J. It is clear that the outboard aileron position is optimum among 
those whose span is 0. 130 of the wing semispan. The improvements over the 
nominal position are less than 5 percent. 

The previous results indicate the aileron should be as far outboard as 
possible. If further improvements are desirable it is necessary to increase 
the aileron span. We held the outboard position at its rra ximum value and 
examined increasing the aileron span by moving the aileron inboard position 
further inboard to the inboard extreme. We did this in two steps -rows 6 
and 7. The values of 6, 6, u, and J are decreased by about 45 percent. 

There are no adverse side effects of using this optimum surface place- 
ment solution on this unflapped vehicle. We established this by looking at 
the roots, gains, and RMS responses. 

Insensitive Design 

Given a dynamic pressure (q) dependent model 
x = A(q) x + BjU + B 2 n 

we want to determine a constant full-state feedback control matrix 
enforces good performance over a large dynamic pressure range, 
the uncertainty weighting procedure (ref. 14). We use a weighting 


(18) 

K that 
We use 
matrix 
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Q = Q + X (6A)' (fiA) 


(19) 


where 

fiA = A(q) - A(q o ) (20) 

q Q = Design dynamic pressure 

q = An off-design dynamic pressure 

X = Control designers insensitivity parameter 

For this application, the design model was obtained by retaining only the 
first two flexure modes, the third-order actuator and the first-order approxi- 
mation to the Dryden wind model. The matrix Q was taken as the case 0 
matrix of Table 1. Design and off-design values of q Q and q are taken both at 
5. 27 and 7. 66 k Pa ; they were also taken at 7. 66 and 5. 27 k Pa. We deter- 
mined that a good value for X is 0. 1. This is the minimum value of X that 
yields about all of insensitivity achievable. 

Baseline and insensitive controllers were designed for two nominal 
values of dynamic pressure q Q = 5. 2 7 k Pa and q Q = 7.66 k Pa. The per- 
formance of the insensitive controllers was compared with that of the nomi- 
nal controllers as the dynamic pressure was varied. Performance evalua- 
tions were based on comparing eigenvalues and RMS responses as dynamic 
pressure was varied between 5. 27 k Pa and 7, 66 k Pa. 

i 

For the designs (baseline and insensitive) developed for 7. 66 k Pa , 
performance was evaluated as dynamic pressure was reduced to 5. 27 k Pa. 
The response produced by the insensitive controller showed a slightly smaller 
variation with dynamic pressure than the baseline controller but the reduction 
in sensitivity was not significant. The baseline controller produced responses 
which did not vary a great deal with dynamic pressure and it was difficult to 
improve upon this performance. 
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For the designs developed for 5. 27 k Pa , performance was evaluated 
as dynamic pressure was increased to 7. 66 k Pa, In this case the baseline 
controller went unstable at approximately 5. 98 k Pa ; whereas, the insensi- 
tive control was stable up to 7, 66 k Pa, Furthermore, the RMS values of 
control surface displacement and velocity did not greatly exceed the allow- 
able values at the higher values of dynamic pressure. In this case the insen 
sitive controller was clearly superior to the baseline controller. 


Control-Observe- Ability 

In the following the quality of controllability and observability are dis- 
cussed pertaining to the control surface/ sensor placement problem. The 
procedure given is due to R. G. Brown (ref. 15), 

Given the nth order system 

x = Ax + Bu (21) 

r = Mx 


The controllability and the observability matrices are given by 


CM 

OM 


A n "^Bl A n " 2 Bl | AB | bJ 

'(a^V |(A T ) n " 2 M T | ~~ | A T M T ] M T ] 


( 22 ) 


Notation: OMCC is CM or OM with columns normalized by normalizing as 

CM or OM is being formed. The magnitudes of the elements of OMCC are 

within the limits of the computer word length and hence produces more 

accuracy. For example, given the k^ stage of computation OMCC of CM, 
s t 

the (k+1) stage is calculated as follows: 
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OMCC k = [o | — I A k-1 B | A k " 2 B | I Bj (23) 

OMCC k+l/2 = f" O | 1 A A k_1 B | A k_1 B 1 A k_2 B | — ! *b] 

OMCC k+1 = [o I --- | A k B I A k_1 B I ! B 


where the overscore bar indicates column normalization. The m column 
vectors of OMCC are denoted as v. and since OMCC is normalized, v [ are 
unit vectors and v^ e R n . Now consider these m unit (Euclidian) n-vectors 
v^ (n £ m) placed in an n by m matrix OMCC 


OMCC = [vj v 2 • • • vj (24) 

We seek n unit n-vectors u. that extremize the cost function 

3 

L = ( V 1 Tu ]) 2+ ( V 2 Tu j) 2 + (25) 


Result 1. If u is any unit vector, L(u) is cost in the direction u. If 
are an orthonormal basis, then 

n 

T< L ^Uj ^ = m (26) 

3 = 1 

The extremizing unit vectors u. form an orthonormal basis and are the 

3 

eigenvectors of matrix OMDD where 

OMDD = OMCC*OMCC T (27) 

Result 2. The associated eigenvalues are the "costs" in the direction 
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( 28 ) 


L. = u. T * OMDD * u. = X. 

3 3 3 3 

X j * 0 (29) 

n 

^ = m (30) 

1 

If L j ( U j) = °‘ ' U j is nearl y orthogonal to the vectors of OMCC. For unit 
vectors u.. such that ju^j si,, the vectors of OMCC span u. with high 
quality. ^ 

To use these results to evaluate the quality of controllability or observ- 
ability we calculate the controllability or observability matrices, normalize 
the columns and calculate the eigenvectors and eigenvalues of the OMDD 
matrix. We found the assessment to be superior to the rank test but defi- 
cient in shading the quality of control-observe-ability. 

The above results can also be used to assess the quality of a prescribed 
accelerometer complement for spanning an open-loop mode acceleration sub- 
space. We computed how well two accelerometers span the subspace of the 
first two acceleration modes, The best results are presented in Table 3; 
accelerometers at node points 16 and 57 are best. In Section V we will also 
show, using another criterion, that this is an effective choice. 
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TABLE 3. - BEST TWO ACCELEROMETERS FOR TWO MODES 


Minimum 

eigenvalue 

Accelerometer 
node points 

0. 999 

16 

57 

0. 910 

16 

45 

0. 718 

16 

30 

0. 694 

28 

57 

0. 611 

28 

45 

0. 526 

16 

18 

0. 462 

16 

55 

0. 439 

28 

30 

0. 300 

16 

43 



SECTION V 

ROBUST KALMAN ESTIMATOR DESIGN AND 
SENSOR OPTIMIZATION 


With full-state feedback the system is unconditionally stable, has good 
gain and phase margins, and has low bandwidth. As predicted by Doyle (ref. 
16) with the ordinary Kalman estimator the system is conditionally stable, 
has poor gain and phase margins, and has high bandwidth. A procedure due 
to Doyle and Stein (ref. 17) is presented to synthesize robust Kalman esti- 
mators. This procedure permits compromising RMS response activity 
against robustness. Here robustness is characterized as gain and phase 
margins. If the sensor complement is minimum phase under control gains, 
full-state stability characteristics can be asymptotically achieved. The re- 
sults of robust Kalman estimator design are presented for three sets of sensor 
complements. Three levels of robustness are used for each of the sensor 
complements and the results compared with that of ordinary Kalman esti- 
mators. 

An optimization procedure, similar to the one employed in Section IV 
for control surface placement, is used to determine sensor locations that 
minimize the cost and RMS activity. Improvements of 15 percent in RMS 
responses are obtained. 

Finally, the sensor position procedure of Stein and Harvey (ref. 18) 
is modified and applied. This yields the same sensor complement previously 
determined by control- observe -ability considerations. 

Robust Kalman Estimator Design 

As described in (ref. 17), the robust estimator design is obtained by 
adding an extra process noise directly to the control input of the plant during 
the estimator design. 
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The equations of motion are given by: 


x = Ax + B- u + B 1 n + B 0 n (31) 

i ic ^ g 

m = Mx + <n (32) 

m 

Here n and n are the nominal gust and measurement noise with 
'g 'm 

E { T1 g (t)T1 g (T) } = W g 6(t " T) (33) 

E { T 1 m (t)T 1 m (T) } = W m 6(t ‘ t) (34) 

and n is the additional process noise added at the control input with 

E {V t)rT c (T) } = qW c 6(t " t) (35) 


and q is the non-negative robustness parameter. 

The robust estimator dynamics is given by: 

x = Ax + B^u + L (Mx - m) (36) 

The robust estimator gain is given by: 

L = - P M'W " 1 (37) 

e m 

The error covariance matrix, P e> is the solution of the Ricatti equation: 

(A + LM)P + P (A + LM) '+ LW V 
e e m 


( 38 ) 



For q = 0, we get the ordinary Kalman estimator design and as q is increased 
the estimator design becomes more and more robust and approaches the 
stability characteristics of the full-state feedback design. 

The control u is given by 

u = Kx (3 9) 

where the gain matrix K is obtained from the baseline case 0 of Table 1 
(with first-order wind model). 

Results for three sets of sensors are presented. In the first set, two 
accelerometers are used; one each at node points 16 and 57; this choice was 
determined by the control-observe -ability considerations in Section IV. For 
the second and third sets, single accelerometers are used at node points 
18 and 28, respectively. 

Ordinary Kalman estimators (q = 0 in equation (8)) were synthesized 
and the performance is summarized in Table 4. Row 1 presents the results 
with full-state feedback control; rows 2, 6, and 10 present the results with 
the three different sensor complements. 

The values of RMS responses are modestly larger when Kalman esti- 
mators are used. The increases in 6 are 29, 32, and 45 percent, respec- 
tively; 5 increases by 13, 20, and 18 percent, respectively. From this 
we see the first (two sensor complement) is best and that the second (node 
point 18) is preferable to the third (node point 28). All of them comfortably 

meet the objectives of 4 deg and 250 deg/sec for 5 and g. 

\ 

The gain and phase margins of full state feedback control are guaranteed 
by Safanov and Athans (ref. 19); -6 dB and +« dB in gain and ±60 deg in phase. 
Row 1 shows that we comfortably exceed these values. The S(full-state feed- 
back curve) of the Nyquist plot of Figure 11 shows the superlative stability 
characteristics. 
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TABLE 4. - ROBUSTNESS SUMMARY 


Control Sensor 


State State 
Kalman 16,57 
Robust 16,57 
Robust 16, 57 
Robust 16,57 
Kalman 1 8 
Robust 18 
Robust 1 8 

Robust 18 
Kalman 2 8 
Robust 28 
Robus t 28 
Robust 28 


0. 0 

0 . 0 

0. 000001 

0. 00001 

0 . 0001 

0 . 0 

0 . 000001 

0 . 00001 

0. 0001 

0. 0 

0 . 000001 
0. 00001 

0 . 0001 


RMS 

5 5 

Deg 

Deg/sec 

1. 804 

164.9 

2. 330 

186.2 

2. 733 

196. 8 

3. 066 

212.0 

3. 534 

240. 2 

2. 383 

197. 9 

4.018 

158. 5 

6.445 

197.4 

11. 534 

243.4 

2. 607 

194. 3 

3. 781 

164. 6 

5. 069 

188. 2 

6. 998 

204. 7 


Gain margin 

dB Hz dB 1 Hz I, Deg [ 


9,90 “ -63.7 

10.28 + 6. 4 21. 87 -42. 5 

9.85 20.0 75.28 -45.4 

9.87 21.5 103.56 -53.3 

- ,. , 9. 90 24. 9 149. 65 -57. 9 

-3.9 10.49 + 3.0 17.1 -43.9 

-4. 5 9. 82 + 8.4 4. 84 -29. 1 

-5. 6 9. 72 + 7. 7 3. 95 -34. 3 

-6. 3 9. 74 + 7. 9 3. 01 -40. 5 

-1.8 61.20 + 3. 7 18.54 -41.5 

-5.6 4.85 +26.3 2.21 -37.9 

-6.8 9.80 +28.9 1.56 -46.2 

-7.5 9. 85 +50.7 0.37 -55.2 


Phase 

mars 

nn 



Bandwidth 

Hz 

Deg .. 

Hz 

Hz 

7. 

05 

+82. 

3 

15. 

25 

23. 2 

8. 

26 

+32. 

4 

13. 

64 

59. 8 

7. 

46 

+66. 

2 

13. 

56 

19.3 

7. 

23 

+72. 

9 

14. 

23 

11. 1 

7. 

16 

+76. 

3 

14. 

64 

21. 8 

8. 

71 

+18. 

3 

1 13. 

15 

44. 9 

7. 

97 

+43. 

9 

12. 

21 

15.3 

7. 

48 

+5 8. 

2 

12. 

62 

16. 6 

7. 

21 

+65 

. 6 

13. 

. 10 

17. 8 

8. 

. 55 

+22 

. 8 

13. 

. 62 

64. 6 

7 

. 81 

+55 

. 6 

12 

. 75 

16. 7 

7 

. 31 

+71 

. 1 

13 

. 67 

18. 5 

7 

. 10 

+76 

. 6 

14 

. 42 

21. 2 












As predicted by the work of Doyle (ref. 16), Table 4 shows that with 
ordinary Kalman estimators the magnitudes of negative stability margins 
are less than 6 dB and phase margin magnitudes are less than 60 deg. The 
K (Kalman) curves of Figures 11-13 demonstrate the poor stability charac- 
teristics. 

Though not directly addressed by theory, it is interesting to note band- 
widths. We define the bandwidth as the highest frequency at which the gain 
is greater than -6 dB. State control has a bandwidth of 23. 2 Hz; good con- 
sidering the open loop flutter frequency of 10. 0 Hz. With ordinary Kalman 
estimators bandwidths are doubled or tripled. 

Sensors 16 and 28 have no zeroes in the right half plane. Thus, the 
first (16, 57) and third (28) sensor complements meet asymptotic robust- 
ness requirements. Sensor 18 has a zero in the right half plane so robust- 
ness capabilities are expected to be more limited. 

For each of the three sensor complements we used values for q of 
0. 000001, 0. 00001, and 0. 0001. Table 4 summarizes the results. 

The RMS responses are generally monotone increasing with q. For the 
first case the increases are modest; even with q = 0. 0001, the values of 5 
and 6 are acceptably small. For the non-minimum phase case (18), the in- 
creases in 6 ace very large and even with q = 0. 000001 are in excess of ob- 
jectives. For the minimum phase single sensor case (2 8) the increases with 
q are more modest. 

Stability margins are also essentially monotone increasing with the 
robustness parameter q. Although negative phase margins are below objec- 
tives; we select q = 0. 000001 as a compromise between RMS response and 
stability objectives. 

Although not assured by theory it is no surprise that bandwidths are 
much smaller with the robust estimators. 
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Figure 13. - Nyquist plot for node 28. 
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The Nyquist plot of Figure 11 for nodes 16, 5 7 demonstrates the effec- 
tiveness of the robust Kalman estimator synthesis procedure: monotonicity 
of margins and transition from Kalman estimator to full-state feedback as 
robustness increases. Low bandwidth is achieved at the lowest robustness 
level. 

Figure 12 shows that with the non-minimum phase sensor, robustness 
is more difficult to achieve. With the minimum phase sensor (28) the results. 
Figure 13, are much like those of Figure 11. 


Sensor Optimization 

With the full- state feedback control law given, sensor positions are 
sought to achieve best system performance when implemented with a robust 
Kalman estimator. The single sensor and two sensor cases are considered. 

Single sensor. -Previous considerations have amply demonstrated that 
the sensor should have minimum phase. We placed sensors at eight positions 
along the rear spar and two on the front spar and computed zeroes. They are 
presented in Table 5. Not surprising, considering ILAF (identical location 
of accelerometer and force), the three node positions (16, 19 and 22) on the 
rear spar and centered about the aileron are minimum phase. However, 
it was not anticipated that the range over which minimum phase is achieved 
would be so small. 

The performance results of full-state feedback control are tabulated in 
row 1 of Table 6. Robust Kalman estimators with q = 0. 000001 are synthe- 
sized with sensor at node points 16, 19 and 22, and the system performance 
evaluated are listed in rows 3-5 of Table 6. Node point 16 is better than node 
point 19 which is better than node point 22. This ordering correlates with 
the sensor zeroes (Table 5) meaning that more negative is better. The gain 
margins are good, but the phase margins need 15 to 20 deg increases in 
magnitude. The control surface deflection 6 exceeds the objectives of 4. 0 
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deg. Rates (5) are below the objective of 250 deg/s. Increasing the robust- 
ness parameter q would increase the phase margins at the expense of RMS 
surface activity; and vice versa. So, we seek to improve performance by 
moving the sensor. We have exhausted the rear spar positions so we posi- 
tion them aft of the spar. The notation 19 (2T) indicates the point twice as 
far from the elastic axis as the node point 19 (Figure 14 shows the DAST 
ARWl node points). 

table 5. -two most positive zeroes* 


Sensor 

node 

position 



2 

+5 92. 0 

+592. 0 

4 

+462. 6 

+462. 6 

7' 

+343. 0 

+343. 0 

10 

+183. 3 

+183. 3 

13 

+ 50. 8 

+ 50. 8 

16 

- 9. 1 

- 9. 1 

19 

- 7. 5 

- 7. 5 

22 

- 6. 1 

- 6. 1 

6 

732. 8 

+971. 2 

12 

679. 8 

+6 7 9. 8 


Neglecting both the two at zero 
and the two wind zeroes at -4. 5. 


The last seven rows of Table 6 present the results. Moving aft and 
toward the tip improves performance. Although 4(3T) has better performance 
than 4(2T), we choose 4(2T) as our final result; 4(3T) is not used because we 
are concerned about accessibility. 

Using the node point accelerometer 4(2T), gain margins are excellent, 
but we would like about 7. 3 deg more phase margins. The surface rate and 
deflections meet the design specifications. 
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TABLE 6. SINGLE SENSOR ROBUST CONTROLLER PERFORMANCE 


Sensor 

Control I node 
position 


RMS per fps gust 


Deg D eg/ sec 


Gain margin Phase margin Bandwidth 


Hz dB Hz Deg [ Hz I Deg I Hz Hz 


State 


0.0 

2.579 

189. 4 

2. 630 

0. 5529x 10 6 

-14. 2 

38. 44 

8. 1 

9. 77 

-65. 1 

6. 72 

95. 1 

14. 21 

86. 3 

Kalman 


4(2T) 

0. 0 

3. 333 

186.6 

3.350 

0. 8565 x 10 6 

- 3. 7 

23. 76 

7. 3 

10.44 

-65.4 

7. 67 

31. 9 

15.64 

28. 0 

Rob us t 


16 

0. 000001 

4. 047 

218.7 

4.584 

0. 1385x 10 7 

-14. 9 

93. 60 

5. 7 

9. 99 

-43. 3 

7. 93 

48. 6 

12. 80 

16.2 

Robust 


19 

0. 000001 

4. 156 

223. 7 

4. 658 

0. 1484 x 10 7 

-15. 1 

94. 33 

5. 5 

9. 99 

-41. 9 

7. 94 

48. 0 

12. 70 

16. 1 

Robust 


22 

0. 000001 

4. 201 

225. 1 

4. 693 

0. 1444 x 10 7 

-15. 1 

95.41 

5.5 

9. 98 

-41. 7 

7. 94 

48.0 

12. 68 

16. 1 

Robust 


1 9(2T) 

0. 000001 

3. 888 

214. 8 

4.508 

0. 1361 x 10 7 

-14. 6 

102. 56 

6. 3 

10. 03 

-50. 1 

7. 72 

52. 8 

13. 05 

16. 6 

Robust 


16(2T) 

0. 000001 

3. 788 

207. 8 

4.391 

0. 1302 X 10 7 

-14. 9 

101. 41 

6. 4 

10. 04 

-51. 7 

7. 70 

52. 7 

13. 13 

16. 9 

Robust 


13 (2T) 

0. 000001 

3. 732 

207. 8 

4.403 

0. 1310 x 10 7 

-14. 5 

101.41 

6. 4 

10.04 

-53. 0 

7. 66 

53.4 

13.16 

16. 9 

Robust 


10(2T) 

0. 000001 

3. 683 

207. 0 

4. 371 

0. 1295 x 10 7 

-14.5 

99.52 

6. 4 

10. 04 

-52. 2 

7. 67 

53. 1 

13. 19 

16. 9 

Robust 


7(2T) 

0. 000001 

3. 630 

207. 9 

4.336 

0. 1283 x 10 7 

-14.5 

98. 01 

6. 5 

10. 06 

-52.4 

7. 67 

52. 8 

13.23 

17. 0 

Robust 


4(2T) 

0. 000001 

3. 595 

209. 7 

4.314 

0. 1277 x 10 7 

-14. 6 

97.43 

6. 5 

10. 06 

-52. 7 

7. 66 

52. 7 

13.24 

17.0 

Robust 


4(3T) 

0. 000001 

3. 471 

205. 4 

4. 232 

0. 1253 x 10 7 

-14. 3 

105. 52 

7. 0 

10. 07 

-57. 3 

7. 37 

56. 2 

13. 48 

17.2 






























Y-AXIS (INCHES) 


X-AXIS (INCHES) 
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A comparison (rows 2 and 11) of the results with the Kalman estimator 
and robust estimator for node points 4(2T) shows robustness cost 8 and 13 
percent in deflection and rate. It improved the negative phase margin 
from the undesirable value of -3. 7 to the desirable -14. 6. Similarly, the 
positive phase margin increased from 31. 9 to 52. 7 deg and the bandwidth 
was reduced from 28. 0 to 17. 2 Hz. Using accelerometer positions 22 and 
4(2T) as extremes, it is worthy of noting the benefits achieved by searching 
for better sensor positions. We improved 5 by 14 percent, 5 by 7 percent 
and J by 12 percent. We also improved by 1. 0 dB in positive gain margin, 

9 deg in negative phase margin and 5 deg in positive phase margin. 

Two sensors. -Here two methods are used to search for optimum sensor 
positions. The first is optimization of sensor location similar to that 
employed for the single sensor case. The second is a modification of the 
Stein-Harvey procedure (ref. 5). 

For robust estimator design with two sensors and a single control, it is 
sufficient if some linear combination of the two sensors is minimum phase. 
Rather than attempt this calculation, we assume that if one of the 
accelerometers is near the aileron, then we can get minimum phase. As 
in the single sensor case, we start from a few promising positions and then 
perturb the positions to determine whether we are near a best position. 

The first 20 rows of Table 7 present our results. The first row 
contains the full state-feedback results. Kalman and robust estimator data 
with sensors at nodes 2 and 6 are presented in rows 2 and 3. We consider 
these to be the best positions for sensors restricted to the front and rear 
spars. The robust aileron RMS value of 3. 008 is considerably less than 
the objective of 4. 0 and the RMS value for rate of 224. 4 is slightly smaller 
than the objective of 250. The gain margins and bandwidth are good. We 
would like to increase the negative and positive phase margins by 20 and 8 
deg, but believe the tradeoff with RMS responses to be undesirable. 
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TABLE 7. TWO SENSOR ROBUST CONTROLLER PERFORMANCE 

05 


Control 

Sensor 

node 

positions 

q 

RMS per fps gust 

J 

Gain margin 

Phase margin 

Bandwidth 

6 

a 

u 

Deg 

Deg/sec 



dB 

Hz 

dB 

Hz 

Deg 

Hz 

Deg 

Hz 

Hz 

State 


0.0 

2.579 

189. 4 

2. 630 

0. 5529 xlO 6 

-14. 2 

38.44 

+8. 1 

9. 77 

-65. 1 

6. 75 

95. 1 

14. 21 

86. 3 

Kalman 

2, 6 

0. 0 

3. 033 

203. 6 

3. 073 

0. 7501x 10 6 

- 2. 7 

21. 44 

+5. 6 

10. 42 

-48. 1 

8. 28 

27. 2 

14. 66 

62. 1 

Robust 

2, 6 

0. 000001 

3. 008 

224. 4 

3. 894 

0. 105 lx 10 7 

-15. 0 

90. 99 

+5. 7 

9. 92 

-40.4 

7. 80 

52. 1 

12. 83 

16. 6 

Robust 

40, 42 

0. 000001 

3. 748 

264. 0 

4. 328 

0. 1325X10 7 

-12. 7 

3. 98 

+5. 1 

9. 85 

-34. 4 

7. 86 

48. 2 

12. 38 

15. 3 

Robust 

31,33 

0. 000001 

3.593 

256. 7 

4. 218 

0. 1238x 10 7 

-15. 0 

3.52 

+5. 3 

9. 85 

-36. 6 

7. 80 

51. 0 

12.53 

15.6 

Robust 

22, 27 

0. 000001 

3.402 

242.4 

4. 082 

0.1151X10 7 

-14. 7 

94.49 

+5. 6 

9. 87 

-38. 9 

7. 75 

53. 2 

12. 68 

15. 9 

Robust 

19, 21 

0. 000001 

3. 430 

234.3 

4. 099 

0. 1150 x 10 7 

-14. 6 

94. 51 

+5. 6 

9. 88 

-39. 0 

7. 78 

52. 6 

12. 68 

15. 9 

Robust 

19, 21 + 18 

0. 000001 

3.446 

231.2 

4. Ill 

0.1151 xlO 7 

-14. 6 

94. 54 

+5. 6 

9. 88 

-39. 0 

7. 78 

52. 5 

12. 68 

15. 9 

Robust 

19, 18 

0. 000001 

3. 463 

228.4 

4. 123 

0. 1155 xlO 7 

-14. 6 

94. 62 

+5. 6 

9. 88 

-39. 0 

7. 78 

52. 4 

12. 68 

15. 9 

Robust 

19 ; 16 , is 

0. 000001 

3. 413 

227. 8 

4. 088 

0. 1138x 10 7 

-14. 6 

94. 22 

+5. 6 

9. 88 

-39. 3 

7. 78 

52. 6 

12. 72 

16. 1 

Robust 

16, 18 

0. 000001 

3. 368 

227. 0 

4. 057 

0. 1124 xlO 7 

-14. 6 

93. 95 

+5. 6 

8. 40 

-39. 6 

7. 78 

52. 8 

12. 73 

16. 2 

Robust 

16, 57 

0. 000001 

3. 232 

237. 5 

3. 988 

0. 1115 XlO 7 

-14. 6 

93. 68 

+5. 8 

9. 88 

-40. 7 

7. 73 

54. 3 

12. 83 

16.4 

Robust 

16, 21 

0. 000001 

3. 345 

232. 3 

4. 041 

7 

0. 1122 x 10 ' 

-14. 6 

93. 85 

+5. 6 

9. 88 

-39. 5 

7. 77 

53. 0 

12. 73 

16. 1 

Robust 

10,12 

0. 000001 

3. 278 

220. 6 

3. 988 

0. 1087 xlO 7 

-14. 8 

91. 93 

+5. 7 

9. 90 

-40. 0 

7. 80 

52. 3 

12. 78 

16. 4 

Robust 

4, 6 

0. 000001 

3. 200 

222. 9 

3. 926 

0. 1062 xlO 7 

-14. 9 

91. 72 

+5. 7 

9. 92 

-40. 4 

7. 80 

52. 3 

12.81 

16. 4 

Robus t 

4, 12 

0. 000001 

3. 177 

223.4 

3. 484 

0. 1053 x 10 7 

-15. 2 

89. 72 

+5. 7 

9. 92 

-40. 0 

7. 81 

51. 7 

12. 81 

16. 6 

Robust 

16,33 

0. 000001 

3. 263 

241.4 

3. 985 

7 

0. 1111 xlO* 

-14. 6 

93. 73 

+5. 7 

9. 88 

-39. 7 

7. 73 

53. 7 

12. 76 

16. 2 

Robust 

10,30 

0. 000001 

3. 216 

232. 6 

3. 938 

0. 1090 x 10 7 

-14. 7 

90.59 

+5. 7 

9. 90 

-39. 7 

7. 77 

52. 7 

12. 78 

16. 4 

Robust 

19, 48 

0. 000001 

3. 289 

243. 5 

3. 989 

0. 1119x 10 7 

-14. 6 

94. 14 

+5. 7 

9. 88 

-39. 9 

7. 72 

54. 2 

12.76 

16. 2 

Robust 

4, 6(2T) 

0. 000001 

2. 984 

214.2 

3. 788 

0. 9904 xlO 6 

-14. 7 

105. 09 

+6. 0 

9. 87 

-44. 3 

7. 61 

57. 7 

12. 88 

16. 4 

Robust 

16, 57* 

0. 000001 

4. 320 

231.3 

4. 7 17 

0. 1447 x 1 0 7 

-15. 8 

23.62 

+5. 0 

9. 96 

-36. 6 

8. 04 

44. 1 

12.49 

15. 6 

Robust 

2, 12* 

0. 000001 

3. 706 

240. 1 

4. 388 

0. 1420 x 10 7 

-16. 4 

87. 87 

+6. 5 

10. 11 

-51. 2 

7. 59 

49. 7 

13.32 

17.5 


Blended to form a single sensor with gains given in Table 8. 









1 

Transfer function of sensor 1 

Transfer function of sensor 2 

The transfer function with numerical values is the unstable flutter mode 
that we are trying to match. In this sense, this is a minimum energy pro- 
cedure (ref. 5). 

To use the procedure, we select sensor positions, compute minimizing 
gains g^ and g ^ and evaluate the least square cost J. The relative values off 

J for different sensor combinations are measures of relative merit. We then 

'■ — 

use the most likely candidates to synthesize a robust estimator to evaluate 
the final result. N 

Table 8 shows results for the least squares part of the procedure. Sen- 
sor node positions are presented in column 1. Nodes 16 and 5 7 and their 
perturbations in the next three combinations were selected on the basis of 
control-observe-ability (Section IV). The last six combinations are the best 
candidates based on achieving minimal controller cost. The costs (column 
4) are not markedly different. The first four have lower costs with the com- 
bination at 16, 57 being minimal. Node combination 2, 12 has minimal 
magnitude for gains g^ and g^\ this might be significant when we evaluate 
system performance using realistic sensor noise. 

We selected 16, 57 and 2, 12 for further consideration. Robust esti- 
mators were synthesized and the systems evaluated. Rows 21 and 22 of 
Table 7 present the results. They are worse than most others. The node 
2, 12 combination is better than that for the node 16, 5 7 combination. In 
this case the magnitude of the sensor noise is an important consideration 
for sensor placement. 



u 

m 2 


u 
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TABLE 8. -SENSOR FREQUENCY RESPONSE MATCHING 


Sensor 

node 

position 

g l 

g 2 

£e 2 

16,57 

-0. 4333 

-0. 2516 

0.3139E + 09 

13,57 

-0. 3948 

-0. 3568 

0. 3239E + 09 

13, 54 

-0. 3877 

-0, 3429 

0. 322 7E + 09 

10,54 

-0. 3483 

-0. 4225 

0. 3327E + 09 

2, 6 

-0. 4988 

+0. 3492 

0. 3891E + 09 

2, 9 

-0. 4424 

+0. 2 977 

0. 4047E + 09 

2, 12 

-0. 3569 

+0. 1950 

0. 4183E + 09 

4, 12 

-0. 4129 

+0. 2428 

0. 4027E + 09 

4,6 

-0. 5448 

+0. 3652 

0. 3723E + 09 

4, 6 (2T) 

-0. 4213 

+0. 2412 

0. 3719E + 09 



SECTION VI 

PRACTICAL FLUTTER CONTROLLERS 


The robust Kalman estimator together with the optimal state feedback 
gains constitute the optimal controller that can be implemented. However, 
the cost of implementation can be significantly reduced by developing lower 
order approximations to the optimal controller. The low-order approxima- 
tions that satisfy the performance requirements are called practical control- 
lers. Two methods of developing practical controllers from the optimal 
controller are residualization and frequency response matching (ref. 4 and 5). 
These methods are used to develop practical flutter controllers for the DAST 
ARW1 wing. Also a practical flutter controller is evaluated at various mach 
numbers and dynamic pressures to check its operation throughout the flight 
envelope. 


Practicalization (Simplification) Procedure 

In simplification, we are attempting to approximate the high order 
optimal controllers with a much lower order approximation that 

• Achieves closed loop asymptotic stability 

• Maintains ample gain and phase margins 

• Maintains low band width 

• Does not markedly increase RMS responses 

We tactily assume that if we approximate the low-frequency character- 
istics with an approximation that ultimately has 40 dB rolloff per decade 
(to satisfy high frequency actuator limitations) then we will meet the objectives. 
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The optimal controller is given by the equations (41) and (42). 


x = |a + B^K + Lc)x - Lm (41) 

u = Kx (42) 

where 

A 

x is the state of the robust Kalman estimator 
K is the optimal state feedback gains 
L is the robust Kalman estimator gains 


The low-order approximations are determined as equations (43) and 

(44). 


A 



■Axj + Bm 


(43) 


u = + Dm 


(44) 


where 

A 

Xj is the state of the low order system. 


Two methods of obtaining A, B, C and D matrices are described in 
the following. 


Residualization , -Here the optimal control system is split into two parts 
as shown by equation (45) and (46). 



A A 

"ll 12 


A 

X 1 

^21 A 22 


A 

X 2 

+ [(LC)j 

(LC) 2 ] 


+ [Mi M2] 


A 


“ — 

X 1 


L 1 




A 



_ X 2_ 


- L 2. 



(45) 
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u 


(46) 



In the first part (low frequency portion) the dynamics is retained. In 
the second part (high frequency portion) only the static terms are retained. 
The static terms are obtained by setting x 2 = 0 and are given by 

h - '[ A 22 + Ma + ( LC ) 2)" 1 I A 21 x 1 - L 2 m l (47) 

Using equation (47) the residualized system is given by 

*1 - { A ll + l B l K )l + < LC) 1 ' A 12 [ A 22 + | B 1 K ) 2 + < LC > 2 ]' 1 A 2lK 

-{ L l - A 12 [ A 22 + | B 1 K ) 2 + < LC > 2 )' 1 L 2 } m (48> 

u m { K 1 “ K 2 [ A 22 + ( B 1 K ) 2 + ^ LC '2) A 2l} x l 

+ K 2[ A 22 + | B l K |2 + < LC) 2]' lL 2 m <49) 

For the practical controllers the rolloff requirement is 40 dB /decade. 

In the examples we worked, the direct transfer matrix is small. So we 
often discarded this term and provided a single rolloff filter at 300. rad / 
sec to meet rolloff requirements. Alternatively, we can retain the direct 
transfer and append another filter at 3000. rad/sec. 

Frequen cy response matchin g .-This procedure (as the name implies) 
involves obtaining a low-order transfer function approximation to the optimal 
controller over a frequency band of interest. 

The transfer function between measurement m and control u (from 
equations 41 and 42) is given by 
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u 


(50) 


= KlsI-A-BjK-Lcl^Vum 
A T(s)m 

For the single control variable case with two measurements 
u = [T l( s) T 2 (s>) 

With one measurement, Tg and m ^ are zero. In frequency response 
matching we approximate each T.(s) with a low-order rational transfer 
function T a ^(s). The numerator and denominator orders are selected and 
we can choose to restrict the T ai (s) to be stable and/or minimum phase. 
Then, the coefficients of T . (s) are calculated to minimize 

a l 

J = w( Uk | I T.( Wk ) - T al | Uk ) | 2 (52) 

where 

u k is a set of points in the frequency band of interest 

W|w k j is a weighting function. 

Practicalization Results 

Here, the optimal flutter controller, obtained with a second-order 
Dryden wind model, is simplified to obtain three practical flutter controllers. 
Two of them use a single sensor aft of wing node point 4. The first approx- 
imates the optimal flutter controller by residualization and the second ap- 
proximates it by frequency response matching; the approximate filters are 
of fourth and fifth orders respectively. The third employs two sensors at 
node points 2 and 6 and is obtained from a seventh-order residualization of 
the optimal controller. All controllers meet the RMS requirements. Sta- 
bility margins are reasonably good. 


m„ 


m r 


(51) 
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Single sensor practicalization results. - Table 9 presents the perfor- 
mance results. | Rows 1, 2, and 3 show the state, Kalman and robust 
results discussed previously. Rows 4 and 5 summarize results with our 
best residualized controller. It is obtained by retaining the first two flex- 
ure nodes in the estimator, by deleting the small direct transfer term (No 
D), and by adding a first-order lag with a 1/300 sec time constant to meet 
rolloff requirements. In row 4 the sensor measurement noise is set to 
zero to permit showing its effect on system performance. A comparison 
of rows 3 and 5 show that residualization has reduced RMS activity. The 
RMS deflection of 3. 41 deg and rate of 184, 3 deg/sec are well below the 
objective of 4 deg and 250 deg/sec respectively. The adverse effect of 
residualization is in the reductions in stability margins. The gain margins 
of -5.0 and +5. 5 are below the objective of ±6. 0 dB. We also suffered a 
13 deg of positive phase margin to reduce its value to about 40 deg. The 
increase in bandwidth from 17. 0 to 25.6 Hz is acceptable. 

Similar results were obtained with frequency response matching using 
a second over third approximation. The RMS values are essentially the 
same, gain margins are good, but we lost another 4 deg in positive phase 
margin. 

Figure 15 presents the Bode plots for the state, Kalman, robust, re- 
sidualized and transfer function approximation cases. The figure shows that 
achieving the desired RMS response at the full- state feedback design level 
force excessive bandwidth on the system. Both practical controllers achieve 
good rolloff after 23. 9 Hz. 

Two sensors practicalization results. -Table 10 presents the performance 
summary. Rows 1, 2, and 3 are for state control and control with Kalman 
and robust estimators. Rows 4 and 5 are for our best controller obtained by 
retaining the first two flexure modes and the gust filter (G) of the estimator 
during residualization; the small direct transfer term was dropped (No. D) 
and the rolloff filter added. The effect of residualization was to yield a small 
but acceptable increase in RMS control surface deflection and a decrease in 
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TABLE 9. SINGLE SENSOR PRACTICAL CONTROLLER PERFORMANCE 


Sensor 

Control node 


RMS | 

6 

5 

u 

Deg 

Deg /sec 

Deg 

2. 58 

189. 4 

2.630 

3. 33 

186. 6 

3. 350 

3.59 

209. 7 

4.314 

3. 16 

180.3 

3. 172 

3.41 

184. 3 

3. 346 

3. 23 

184. 3 

3. 246 

3. 42 

189. 3 

3.434 


Gain margin 


dB Hz dB Hz 


™ • Band- 

Phase margin wLdth 


Deg Hz Deg Hz Hz 


95. 1 14. 2 86. 3 

31. 9 15. 6 28.0 

-52. 7 7. 7 52.7 13.2 17.0 

-65.5 8. 1 '39. 6 16. 1 25.6 

-65. 5 8. 1 39.6 16. 1 25.6 

35. 4 12. 9 27.1 

35. 4 12. 9 27.1 


Retained Added 
Components Filters 


State* 

Kalman 

Robust 

RESID* 

RESID 

TRNFN* 

TRNFN 


* Without sensor measurement noLse 


-14.2 38.4 8. 1 9.8 

- 3. 7 23. 8 7. 3 10. 4 

-14.6 97. 4 6.5 10.1 

- 5. 0 23. 7 5, 5 10.2 

- 5.0 23.7 5.5 10.2 

- 8.3 30. 4 6.5 10.4 

- 8.3 30. 4 6. 5 10.4 


1,2, No D 
1,2, No D 
2/3 
2/3 


300 


S+300 


S+300 

300 


S+300 










































TABLE 10. TWO SENSOR PRACTICA L CONTROLLER PERFORMANCE 




rate. Phase margins are improved and are acceptable. Most serious 
is the loss in positive gain margin. 

Figure 16 compares the frequency responses of state, Kalman robust 
and residualized controllers. We see that the Kalman control attempts to 
maintain the high bandwidth of the state control. The bandwidth of the 
robust is much smaller than that for the Kalman, The residualized practical 
controller is quite similar to the robust. 

Practical flutter controllers. - The three final design controllers are 
presented in Table 11. The two single sensor controllers look good. They 
are stable and minimum phase. Poles and zeros are neither excessively 
large nor small. The two sensor controller (obtained by residualization) 
is non-minimum phase. The transfer from sensor 6 has a zero in the 
right half plan. Otherwise, the controller is satisfactory. 

Performance Analysis of the Final Flutter Controller 

The residualized practical controller with the sensor (accelerometer) 
located at a point on the wing twice as far from the elastic axis as node 
point 4 is chosen to be the final practical controller for further evaluation. 
The sensor location is marked on Figure 14 which shows the DAST ARWl 
node points. Its (x,y) coordinates are 2. 022 m and 1,780 m. The flutter 
controller transfer function is given by the equation 

6 a c _ 1. 1961 (s 2 + 116. 7276 s + 7819. 8649) 

K (s + 300) (s 2 + 73. 6256 s + 2829. 1761) 

(s + 22. 91) 

s 2 + 28. 0845 s + 362. 1409) (53) 

where 6 is the input to the actuator and h is the vertical displacement of 

ctC 

the wing at the accelerometer location. The flutter controller was designed 
58 
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TABLE 11. - PRACTICAL CONTROLLER TRANSFER FUNCTIONS 


Method 

Sensor 

Transfer functions 

RESID 

4(2 T) 

(300 ) |s 2 + 2. * 0.6600 * 88.43 S + 88.43 2 ! (S + 22. 91) 

(S+300 ) (S 2 + 2. * 0.6921* 53.19 * S + 53. 19 2 ) |s 2 + 2.* 0.7379 *19. 03 * S + 19.03 ) 

trnfn 

4(2T) 

006855) < 300 > Is 2 + 2. * 0. 1496 * 92.57 *S + 92.57 2 j 

(S + 300) (s 2 + 2. * 0. 4247 * 70. 04 * S + 70. 04 2 ) (S + 17. 65) 


2 

_ . (300 ) I S 2 + 2 . * 0 . 6258 *18.57 *S + 18.57 2 )(s 2 + 2. * 0.9991* 4. 477* S + 4.477 2 )(S + 158.) 

(U.UUSJO/) (s + 3Qo) ^2 +2 _ # 0.6818 *26.73 *S + 26. 73 2 )(S 2 + 2. * 0. 9938 * 7. 498 * S + 498 2 ) (S + 154.9) (S + 3. 973) 


(S + 300 ) IS 2 + 2. * 0.6818 * 26.73 *S + 26. 73 2 ) |s 2 + 2. 4= 0.9938 *7.498 *S + 7 . 498 2 ) (S + 154. 9) (S + 3.973) 
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by considering the DAST half- wing model with one aerodynamic lag term 
(L=l) and 5 modes. In the following analysis the DAST half- wing will be 
represented by different models, (by changing the number of lag terms and 
modes considered) and the flutter controller will be evaluated. In addition, 
the performance of the flutter controller as the mach number and dynamic 
pressure change will be studied. 


Performance with different models . -Table 12 shows the evaluation with 
four different models for DAST half wing. The first model is used during 
the design. We notice that the RMS values go up as more number of modes 
and aerodynamic lag terms are included. The gain and phase margins do 
not change considerably. The bandwidth is about constant except in the case 
of fourth model. The fourth model with number of lag terms equal to four 
and number of modes equal to eight is considered to be the most accurate 
representation of the DAST half wing. The fifth model is identical to the 
fourth except that there is no sensor noise. So the gain margin, phase mar- 
gin, and bandwidth data are identical to that of the fourth model. The differ- 
ence in the RMS values for the control surface between the fourth and fifth 
model is due only to the sensor noise. 


Performance at different mach numbers and dynamic pressures. -Here 
the model used contains four aerodynamic lag terms and five modes. The 
results are presented in Table 13. We notice that the RMS values at any 
Mach number increases as the dynamic pressure is increased. Also, the 
bandwidth of the system increases as the dynamic pressure is increased. In 
terms of stability, the worst gain margin is 2. 91 dB and the worst phase mar- 
gin is -14. 81° (note that both of these are at dynamic pressure of 9. 58 k Pa ), 
If we restrict the dynamic pressure to be less than 7. 66 k Pa , the worst gain 
margin is 3. 03 dB and the worst phase margin is 25. 6°. Table 14 presents 
the fluter dynamic pressure and the flutter frequency for the three Mach num- 
bers. At Mach number 0. 9 the dynamic pressure at which the DAST half 
wing flutters is more than doubled with the flutter controller. 
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TABLE 12. PERFORMANCE VARIATION WITH DIFFERENT APPROXIMATIONS 
TO THE DAST HALF WING (MACH NO. = 0. 9 q = 7. 66 k PAS) 


No. 

Model used 
for DAST 
half wing 

Control surface 
deflection (RMS) 
(Deg) 

Control surface 
velocity (RMS) 
(Deg/sec) 

Gain margin 

Phase margin 

Bandwidth 

dB 

Hz 

dB 

Hz 

Deg 

Hz 

Deg 

Hz 

Hz 

1 

L = 1 ; 5 modes 

3.41 

184.33 

-5. 0 

23. 7 

5. 5 

10. 20 

-65. 5 

8. 10 

39. 6 

16. 1 

25. 6 

2 

L = 1 ; 8 modes 

3.41 

183. 25 

-4. 9 

23. 7 

5. 6 

10. 17 

-63. 1 

8. 12 

39. 1 

16. 1 

25. 6 

3 

L = 4 ; 5 modes 

3. 85 

207. 67 

-5. 5 

23. 8 

5. 21 

9.84 

-55. 36 

8. 11 

36. 7 

16. 2 

25. 6 

4 

L = 4 ; 8 modes 

3. 87 

207.92 

-5.4 

23. 8 

5. 23 

9. 81 

-52. 6 

8. 14 

36. 3 

16. 3 

27. 1 

5 

L = 4 ; 8 modes 

3. 65 

204. 17 

-5.4 

23. 8 

5. 23 

9. 81 

-52. 6 

8. 14 

36. 3 

16. 3 

27. 1 


(No sensor noise) 
























Mach 

Dynamic 

Control surface 

Control surface 


Gain margin 



Phase margin 


Bandwidth 

number 

pressure 
q (k Pa ) 

deflection (RMS) 
(Deg) 

velocity (RMS) 
(deg/sec) 

dB 

Hz 

dB 

Hz 

Deg 

Hz 

Deg 

Hz 

Hz 


3. 83 

1. 79 

84. 63 

-13.47 

23. 83 

- 

- 

-104. 5 

4. 19 

27. 27 

9. 36 

17. 7 


5. 75 

2. 88 

173. 38 

-10. 07 

23. 83 

- 

- 

-112. 76 

4. 76 

27. 12 

10. 57 

18. 5 

0. 7 

7. 66 

4. 44 

284. 88 

-7. 68 

23. 81 

3. 03 

10. 92 

-103. 5 

8. 67 

34.32 

14. 30 

20. 2 


9. 58 

4. 80 

279. 54 

-5.83 

23.77 

2. 91 

10. 84 

-23. 84 

9. 02 

22. 39 

14. 98 

24. 5 


3. 83 

1. 88 

90. 18 

-12. 54 

23.81 

- 

- 

-104. 62 

4. 25 

25. 60 

9. 34 

18. 3 


5. 75 

3. 02 

177. 96 

-9. 15 

23.84 

19. 19 

8. 88 

-113.89 

5.02 

29. 54 

10.38 

19.4 

0. 8 

7. 66 

4. 05 

240. 81 

-6. 78 

23. 84 

3. 95 

10.42 

-73.80 

8. 64 

37. 08 

15. 01 

21. 6 


9. 58 

4. 95 

257. 55 

-4. 95 

23.83 

4. 16 

10. 29 

-22.38 

8. 09 

25.48 

15.93 

25. 6 


3. 83 

2. 03 

97. 93 

-11. 17 

23. 62 

- 

- 

-105. 38 

4. 37 

24. 32 

9. 27 

18. 8 


5. 75 

3. 09 

173. 15 

-7. 85 

23. 77 

13.41 

9. 03 

-114. 27 

5. 48 

37. 88 

10. 15 

20. 7 

0.9 

7. 66 

3. 85 

207. 70 

-5.48 

23. 82 

5. 21 

9. 84 

-55. 36 

8. 11 

36. 70 

16. 23 

25. 6 


9. 58 

6. 07 

274. 91 

-3. 65 

23. 84 

5. 68 

9. 69 

-14. 81 

7. 06 

25.84 

17.44 

27. 1 

















Figures 17 and 18 present the root locus at Mach 0. 9 of the eight modes 
as the dynamic pressure is varied (here the model contains four aerodynamic 
lag terms and eight modes), with and without the flutter controller, 

TABLE 14.- FLUTTER DYNAMIC PRESSURE AND FLUTTER FREQUENCY 
VARIATION WITH MACH NUMBER FOR THE DA ST HALF WING 


Mach 

Without flutter suppression control 

With flutter suppression control 

number 

q f (k Pa ) 

f f (Hz) 

q f (k Pa ) 

f f (Hz) 

0.7 

6. 32 

9. 266 

11. 72 

6. 871 

0.8 

5. 75 

8. 650 

11. 08 

6. 584 

0. 9 

4.96 

8, 029 

10, 39 

6. 272 


The performance analysis indicates that the flutter controller more than 
achieves the required 44 percent increase in flutter dynamic pressure. At 
the design point of Mach number 0, 9 and dynamic pressure of 7. 66 k Pa , the 
RMS activity of the control surface deflectipn and rate are within the design 
objectives. Also the gain margin, phase margin and bandwidth are satis- 
factory. However, over the range of mach numbers and dynamic pressures 
considered, the gain margins and phase margins are not entirely satisfactory 
(with constant gain flutter controller). This indicates the need for gain sched- 
uling of the flutter controller over the range of operation. 
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Figure 18.- Root locus at M = 0. 9 (flutter suppression control). 



SECTION VII 

CONCLUSIONS AND RECOMMENDATIONS 


The objectives of this study were four fold: 

1 Develop an active flutter control methodology 

2 Design a flutter controller for DAST ARW1 wing to increase 
the flutter speed by 20 percent 

3 Investigate methods for surface/ sensor placement, and 

4 Update NASA owned KONPACT software with respect to model- 
ing procedures and design methods developed during the course 
of this study. 

These objectives were primarily met. The active flutter control metho- 
dology, the summary of design results for the DAST ARW1 wing and the 
results of surface/ sensor placement are presented in this volume. The de- 
tailed technical results, documented as customer engineering letters during 
this study, have been collected as appendices and presented in Volume II. 

The modifications to KONPACT software are documented in Volume III. 

In the following, the results and recommendations for future work in 
the area of active control technology are presented. 


Significant Results 

• An active control methodology was successfully applied to design 
a flutter suppression control for DAST ARW1 wing to increase the 
flutter speed by 20 percent. 
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• Effective procedures were investigated and applied to the surface/ 
sensor placement to improve the performances significantly. 

• An automated procedure was implemented to model the unsteady 
aerodynamics of the DAST ARW1 wing and effectively interface 
KONPACT with the NASTRAN software. 

Recommendations for Future Work 

• Develop and apply an integrated active control methodology to 
meet the multiple design criteria which include flutter suppression, 
rigid body control specifications, and other demands (such as 

maneuver load control and gust load alleviation). 

/ 

• Implement algorithms for application of symmetric singular 
value analysis (ref. 19) to study the robustness of multiloop 
systems. 

• Implement algorithms to develop low order transfer functions 
(having the same denominator) to match the frequency response of 
multi input multi output system. This will simplify the task of 
obtaining minimal order practical controllers. 

X 

• Develop efficient algorithms to solve the constrained optimal 
control problems to enforce conventional design criteria into 
optimal active control synthesis techniques. 

• Develop design procedures using residualized Riccati equations 
via singular perturbation and boundary layer techniques (ref, 

20, 21) to reduce design cost for high order systems. 

• Add optimal digital control analysis and design algorithms to 
KONPACT software. 
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• The existing optimal control routines (particularly FFOC: Fixed 
Form Optimal Control) should be revised for higher computational 
efficiency. 

Conclusions 

• Active control methodology, based on linear quadratic Gaussian 
theory, is viable for design of flutter suppression systems for 
flexible vehicles. 

• Optimization of surface /sensor placement offers significant per- 
formance improvement. 

• Use of a lower order approximate model for the flexible wing 
during the flutter control design and verifying the performance of 
the final controller with a higher order more exact model of the 
flexible wing reduces the computer cost of design significantly. 

• Present design objective, for the DAST ARW1 wing, relative to 
RMS and stability requirements are too demanding. More aileron 
effectiveness (by increasing aileron span) would provide the flutter 
control designer freedom necessary to allow design of a practical 
flutter controller that meets the RMS and stability requirements 
more easily. 
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